NASATMX-71810 


NASA TECHN ICAL 

MEMORANDUM 


i 


o 

I— I 

oo 


X 



< 

oo 

< 


(N ASA-TM-X-7 1 8 1 0) ANALYSIS AND CORRECTION 
OF GROUND REFLECTION EFFECTS IN MEASURED 
NARROWBAND SOUND SPECTRA CEPSTRAL 

TECHNIQUES (NASA) 17 p HC $3.25 CSCL 20A 




N76-10820i 


Unclas 

03055. 


ANALYSIS AND CORRECTION OF GROUND REFLECTION EFFECTS IN 
MEASURED NARROWBAND SOUND SPECTRA 
USING CEPSTRAL TECHNIQUES 

by J. H. Miles, G. H. Stevens, and G. G. Leininger 
Lewis Research Center 
Cleveland, Ohio 44135 

TECHNICAL PAPER to be presented at Ninetieth Meeting of the 
Acoustical Society' of America 

San Francisco, California, November 4-7, 1975 ' . 

Reproduced by ' ' . ' 

NATIONAL TECHNICAL 

INFORMATION SERVICE ' > 

US Department of Commerce - . _ . ^ 

, Springfiefd, VA. 22151 _ , > 




5 



E-8503 


ANALYSIS AND CORRECTION OF GROUND REFLECTION EFFECTS IN MEASURED 
NARROWBAND SOUND SPECTRA USING CEPSTRAL TECHNIQUES 
by J, H. Miles, G. H» Stevens, and G. G. Leininger' 

ABSTRACT 


Ground reflections generate undesirable effects on acoustic measure- 
ments such as those conducted outdoors for jet noise research, aircraft 
certification, and motor vehicle regulation. This paper shows how 
cepstral techniques developed in speech processing can be adapted to 
identify the echo delay time and to correct for ground reflection ef- 
fects. A sample result is presented using an actual narrowband sound 
pressure level spectrum. The technique can readily be adapted to exist- 
ing fast Fourier transform type spectrum measurement instrumentation to 
provide field measurement of echo time delays. 
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integer constant 

direct signal power spectral density 

amplitude of the complex reflected signal transfer function 
quefrency 

sound pressure level re 2xl0~^ Pa, dB 
reflectance 
constant 

ratio of path length of the reflected signal to the path length 
of the direct signal 
phase angle parameter, ojx - 9(a)) 
delta function 

phase angle of the complex reflected transfer function 
time delay 

angular frequency, 2irf 
Subscripts; 
m measured 

2 intensity - summed 

Superscript : 

+ without echo effects 

INTRODUCTION 

The fundamental problem discussed herein is to identify and remove 
ground reflection effects so that the true narrowband spectral emission 
properties of a noise source may be determined. The particular problem 
discussed commonly occurs when noise measurements of predominantly low 
frequency noise sources such as full scale jets, large-scale blown flaps 
or combustors are made over a ground plane. Due to NASA activity in air- 
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craft noise reduction, considerable effort at NASA Lewis Research Center 

has gone into experimental and theoretical investigations aimed at mini- 

1-4 

mizing reflection effects. Ground reflection effects are also impor- 
tant because they may cause spectrum measurement errors that can have a 
large economic impact due to various federal, state, and city noise con- 
trol regulations.^’^ 

An analytical solution to the problem for plane waves is simple in 
principle. However, real effects such as distributed sources, nonplanar 
waves , temperature and wind-created shadow zones complicate the problem 
and render analytical solutions difficult. Thus, emphasis has been 
placed on experimental procedures to eliminate ground reflection effects 
such as using ground microphones or overhead-boom microphones (which can 
.be subject to operational problems), and computer model-inverse filter 
type procedures to correct ground reflection effects in measured data 
(which are partially empirical) . 

A new procedure to identify and remove reflection effects is pre- . 

sented herein which avoids previous problems. The technique described 

draws upon an analytical method of signal processing known as cepstrum 

analysis. In the present work, the cepstrum is the Fourier transform of 

the sound pressure level spectrum which has the necessary property of 

being a logarithmic function of the sound pressure. The cepstrum tech- 

7-12 

nique has been used for echo detection, speech analysis (i.eo, find- 

13 14 

ing the pitch period of a voiced speech segment), * and photographic 
15 16 

processing. * This paper discusses how the cepstrum technique was 
adapted to identifying the echo time delays on the ground reflection 
problem. The method is sufficiently general to eliminate phase cancel- 
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lation and reinforcement effects due to multiple reflecting surfaces and 
produces a spectrum resulting from only intensity addition of the direct 
and reflected signals. 

This paper contains first a discussion of the cepstral technique as 
applied to the analysis of ground reflection effected sound pressure 
level spectral datao The calculation method is next discussed, and an 
example is presented using an actual experimental sound pressure level 
spectrum. Finally, some remarks on possible application to real time 
data analysis are given. 

THEORY 

Effect of Reflections on the Spectrum 

The procedure described herein for identifying and removing reflec- 
tion effects is based on the effect reflections have on the spectrum^' 

For a signal having a deterministic waveform or a signal from 'a real ran- 
dom ergodic process, the effect is the same. In Ref. 3, this effect is 
derived for a multiple noise source treated as a real random ergodih 
process® The resulting equation shows that the power spectral density-of 
the received signal H(oj) is the product of the spectral density of the 
direct signal P((a)) and a transfer function corresponding to the reflec- 
tion process called the reflectance |t(o))|^. Thus, 

H(W) = (1) 

‘The reflectance for a single point source is given by 

ItCoj) 1^ = 1+ -2-^ + cos[o)T - 0(o))] (2) 

z 

where y) is the angular frequency,, Q((jj) is the amplitude, and 6(a)) is 
the phase angle of the complex reflected signal transfer function (be- 
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cause amplification does not occur along the reflected signal path Q(o)) 
has a maximum value of 1), t is the echo time delay from source to micro- 
phone and z is the ratio of the path length of the reflected signal to 
the path length of the direct signal. 

The reflectance for a point source (Eqo 2) and the more general re- 
flectance derived in Refo 3 for a distributed source. have a cosinusoidal 
form. This functional form is the basis of the procedure developed to 
correct the measured spectrum. The procedure will be discussed next. 

The Cepstrum . 

As indicated, the power spectral density of the composite signal 
(consisting of the direct and reflected signal) is the product of two 
terms. The first term (P(tjo)) is the power spectral density of the direct 
signal. . The second term (|t(co)| ), due to the echo, is periodic and pro- 
duces a. nearly cosinusoidal ripple in the spectrum of the composite sig- 
nal. This ripple is readily observable in the sound pressure level spec- 
trum which is defined as ten times the logarithm of the power spectral 
density of the received signal. From Eq. 1 the received ;spectrum is 
written as 


SPL(u) = 10 log^QH(o)) = 10 log^QP(m) + 10 logj^QlT(w)r ‘ (3) 

Substitution of Eq. 2 into Eq. 3 produces the expression 

SPL(to) = 10 log^QP(aj) + 10 logj^Q^l + cos[o)T - e(aj)]| (4) 

The effect of reflection on the sound pressure level spectrum is 


periodic. From Eq. 4, the reflectance term has spectral peaks at fre- 
quencies which are integral multiples of 1 /t plus a phase factor 
0(u))/2irTo 
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The ripple observed in the spectrum received by a microphone near 

the ground is from Eq. 4 periodic in frequency. For functions periodic 

in time, the reciprocal of the period (time spacing between cycles) is 

frequency and has dimensions of cycles per second. By convention, to 

prevent any confusion with this latter terminology, the reciprocal of the 

13 14 

period of the spectral ripple is called the quefrency * and has the 
dimension of time. Thus, short period ripple in a sound spectrum corre- 
sponds to high quefrency components. The observed spectrum received by a 
microphone near the ground generally consists then of high quefrency com- 
ponents caused by reflections superimposed on the direct sound spectrum. 

Just as the Fourier transform of a signal can be used to identify 
frequency components of a signal, the Fourier transform of the sound pres- 
sure level spectrum can be used to identify quefrency components in the 
spectrum. Such an analysis produces the cepstrum. While frequency is 
the , independent variable in the spectrum, quefrency, q, is the independent 
variable in the cepstrum. For noise data measured near a ground plane, 
the cepstrum contains large peaks at the quefrencies related to echo de- 
lay times and multiples thereof. Thus, cepstral analysis is an effective 
way to find echo delay times. 

This c^ be shown as follows. First, Eq, 4 is written as 


SPL(co) = SPL' 


,■^(03) + 10 log^Q 1 + 

L z . 




+ 10 log^Q^l + cos[tJT - e(o))]) (5) 




1 + ^ 
2 

z . 


where 
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SPL'*'(a)) = 10 log^QP(oj) 


( 6 ) 


Now, let 



and 

a = toT - 6 (to) (8) 

From the definition of Q(to) and z, it follows that a^ is less than 
or equal to one. Thus, the third term of Eq. 5 can be expanded by the 
series expression for the Naperian logarithm for (1 + y) 


ln(l + y) 



/ 1 \ n— 1 n / 

(-1) y /n 


(9) 


which converges provided 
-I < y <1 

Substitution of Eqs. 7 and 8 into Eq. 5 and using Eq, 9 and the rela- 
tion of logarithms of base 10 with logarithms of base e permits the 
third term of Eq. 5 to be expressed as follows: 


10 log^gCl + a^ cos a) = 10 log^^ge 


n=0 


n 

cos g 
n 


( 10 ) 


Powers of cos a can be expressed as 
(^2 a ^ ^-2.a^ 

cos c 


n _ + e"^ 

’ “ \ 2 ) 


( 11 ) 


Using this relation in Eq. 10 and substituting the result in Eq. 5 shows 


that 
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oo 

SPL(uj) SPL^(w) + 10 log^Q|l + ^ ^ ^ cos^[ioT - 0(w)]^ (12) 

n=0 


The Fourier transform of Eq. 12 is formed to obtain the cepstrum 
C(q). For the case where 0(a)) is a constant times o) and a^ is inde- 
pendent of 0 ), 

C(q) = F[SPL(o))] = F[SPL‘^(o))] 

+ f[i0 

Therefore the cepstrum for the single echo case consists of the cepstrum 
of the "intensity-summed" spectrum plus a train of impulses (i.e., delta 
fi:^nctions) occurring at times which are integral multiples of x - 0(a)) /o). 
Intensity summing implies that the combined energy of the direct and re- 
flected signal is obtained and that reinforcements and cancellations due 
to signal phasing are absent. The intensity-summed cepstrum i-s 

therefore 


C^(q) = F[SPL'^(a))] + F 



(14) 


Intensity-Summed Spectrum Extraction Procedure 
The intensity-summed spectrum is extracted from the cepstrum by 
first smoothing the cepstrum. This is done to remove the delta functions 
due to the echoes and may be accomplished by interpolating between points 
adjacent to the peaks which appear in the cepstrum due to' the echo. The 
resulting intensity-summed cepstrum is then inverse-Fourier-trans formed. 
From Eq. 14 the resulting intensity-summed sound pressure level spectrum 
is 
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SPlJ( 03) - F"^[cJ(q)] = SPL'^Co)) + 10 j (15) 

Thus, all phase cancellations and reinforcements due to echoes have been 
removed, . 

COMPUTATION METHOD 

This section describes the computational procedure used. The f ol- , 
lowing section discusses results obtained with real datao 

Fourier Transform Computation 

The Fourier transform used to compute the cepstrum was performed 
using the fast Fourier transform method of computing the discrete Fourier 
transform. (Fast Fourier transform methods of analyzing sampled time 
series data are described in Refs. 17 to 23.) The Fortran computer pro- 
gram used herein to implement the fast Fourier transform is described in 

k 

Ref. 20 and requires 2 points, where k is some integer. 

Cepstrum Smoothing Procedures 

Methods for smoothing the cepstrum are discussed in Ref. 8, The 
basic procedure is designed to remove the periodic ripple in the spectrum 
introduced by the multiple echoes. This can be achieved by multiplying 
the cepstrum by an appropriate function. The method used herein is 
called ”comb filtering." 

Comb filtering is performed by multiplying the cepstrum by a func- 
tion that is unity everywhere except at those points which are considered 
to be due to echoes. At these points the filter function is zero. The 
resulting filter cepstrum is then found at* the zero points by linearly 
interpolating between the preceding and subsequent nonzero points in the 


cepstrum. 
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Selection of the cepstrum quefrencies due to echoes must be done 
with caution since not all nonzero values of the cepstrum are due neces- 
sarily to echo time delays. For example, the cepstrum of physically 
separated but coherent sources can resemble one with echoes. Also, a 
double peaked spectrum (typical of a suppressor nozzle) will introduce 
nonzero components in the cepstrum which as previously discussed are sub- 
ject to misinterpretation. The proper quefrency selection can be achieved 
in a number of ways. , With the above cautions in mind the simplest proce- 
dure is to assume that all spikes larger than a given value. at quefrencies 
larger than the lowest reasonable echo delay time are actually due to, 
echoes. 

The steps fpllowed. by the method are as follows; 

(1) Form the new sound pressure level spectrum array; : ; 

(2) Calculate the transform to produce C(q); 

(3) Apply the comb filter to produce C^(q) 

(4) Inverse transform the cepstrum to produce SPL^(ai) ; . 

This is schematically diagrammed in Fig, 1. 

. RESULTS AND DISCUSSION . . • 

Application to Measured Data 

Spectrum . - To test the cepstrum procedure on real measured data a 
set, of graphical narrowband jet noise data was selected,. ^ The data are. 
subsonic hot jet noise data from a 10.16 cm nozzle operated with a veloc- 
ity of 548.64 m/sec at a tenq)erature of 649° C (1200°. F) taken over an 
asphalt surface with geometry as follows: -microphone height, 1.676 m; 
source-microphone distance, 7.620 m; source height, 1.67 m; ratio of echo 
path length to direct signal path length, 1.09252 m. One hundred points 
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were read from the narrowband plot. For the cepstrum procedure, 

128 data points (k = 7) were needed and they were obtained by linearly 
interpolating 128 points from the 100 measured points. The first 100 
interpolated points are shown in Fig. 2. The ripple due to ground re- 
flection effects is very pronounced below 4700 Hzo Beyond 4700 Hz the 
ripple effect almost disappears. 

Cepstrum . - The cepstrum was calculated using the method described 
above. Rather than showing the resulting cepstrum which has an approxi- 
mate 4 decade range, the cepstrum level given by 

= 10 loglc(q)!^ (17) 

is presented in Fig. 3. Most of the energy, as discussed earlier, is at 

a quefrency of. zero which corresponds to the energy , associated with the 

direct signals However, significant energy is also apparent at que- 

-2 -2 

frencies between 0.2x10 and 0.205x10 seconds and. again. between . 

-2 -2 

0.405x10 and 0.410x10 seconds. Thus the echo time delay correspond- 

-2 

ing to the large peak at a quefrency of 0.205x10 seconds is easily 
identified by the method. 

Intensity-summed spectrum . - The comb filter was applied to the . 
cepstrum. The points used in the comb filter are indicated by the tick 
marks in Fig. 3. .. 

The resulting intensity-summed spectrum obtained by inverse trans- 
forming the comb filtered cepstrum is shown by the solid curve in Fig. 4. 
The ripple at low frequencies has been> essentially removed. 

Application to Real Time Data Analysis 
The cepstrum method of analyzing sound pressure level spectra, de- 



scribed herein can be used to identify echo time delays and to remove 
reflection effects in measured spectra as the data are being taken. The 
example given herein used graphical data for convenience. Digital sig- 
nal processing instrumentation using fast Fourier transform techniques 
could be modified to produce cepstrum level plots so that echo time de- 
lays could be identified in real time. The comb filtering process and 
inverse Fourier transform calculation could also be performed automati- 
cally so that intensity summed sound pressure level spectra could be cal- 
cula.ted in real time. 

CONCLUDING REMARKS 

This paper has shown that the effect of ground reflections on 
narrowband jet noise spectrum data can be adequately removed by cepstral 
t;echniques. Cepstrum calculation is a signal processing method developed 
originally for speech analysis and used for echo detection and photographic 
processing. The method developed herein extends the use of cepstral tech- 
niques into the field of aeroacoustics. 

When the method was applied to measured jet noise data, it was found 
to eliminate all significant phase cancellation and reinforcement effects 
due to reflecting surfaces, producing a spectrum due to intensity addi- 
tion only of the direct and reflected signals. 

The method could be adapted to current instrumentation to provide a 
direct measurement of time delays due to reflection effects and thereby 
identify sources of echoes. This would be useful for aircraft noise 
certification, motor vehicle noise regulation, arid jet noise research 


acoustic testing. 
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Figure 1. - Diagram of cepstrum 
analysis. 



FREQUENCY, f. KHz. 

Figure 2. - Measured sound pressure level spectrum of jet 
noise. 



Figure 3. - Cepstrum level calculated from spectrum of 
figure 2. 


CEPSTRAL CALCULATED INTENSITY 



FREQUENCY, f. kHz 

Figure 4. - Comparison of measured spectrum and cepstral 
calculated intensity summed spectrum. 
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